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Abstract - Discrepancies of range between the Space Surveillance Network radars and the Deep 
Space Network in tracking the 1998 earth flyby of NEAR, and between ESA’s Doppler and range 
data in Rosetta’s 2009 flyby, reveal a consistent excess delay, or lag, equal to instantaneous one-way 
travel time in the telemetry signals. These lags readily explain all details of the flyby anomaly, and 
are shown to be symptoms of chirp d’Alembertian travelling wave solutions, relating to traditional 
sinusoidal waves by a rotation of the spectral decomposition due to the clock acceleration caused 
by the Doppler rates during the flybys. The lags thus relate to special relativity, but yield distance 
proportional shifts like those of cosmology at short range. 


The fourth-power power law limits direct radar tracking, 
as provided by the Space Surveillance Network (SSN), to 
about the range of geostationary orbits (36,000 km). For 
tracking spacecraft on deep space missions, NASA’s Deep 
Space Network (DSN) uses the telemetry signal returned 
by the phase-coherent onboard transponder for both range 
and Doppler measurements, using modulated range codes 
and the carrier, respectively, as detailed in in §ni] ■ Using 
spin-stabilized spacecraft, this approach achieves sufficient 
precision for tests of general relativity EH!]- Over decades, 
this approach has led to four space anomalies [5], of which 
the best known, the Pioneer anomaly, has now been traced 
to an overlooked radiation reaction [6]. 

The present work fully explains the earth flyby anomaly, 
without assuming dark matter (cf. ,7]), or modifications 
to gravitation theory (cf. [5]f5] ). A broader result is a local 
mechanism that relates more closely to special relativity 
and propagation, yet yields distance proportional spectral 
shifts along with time dilations, which are thought to need 
an expanding space-time (cf. |T0Hl3l ). 

The distance proportionality is given by large negative 
residuals of the SSN data [14] , against the DSN-estimated 
trajectory, which, barring contrived hypotheses, can only 
mean either that the SSN radar echoes were superluminal 
specially during the flybys, or that the DSN Doppler and 
range data had an excess delay. These residuals have been 
omitted in later discussions [9 HT5H23] , as they exceed the 
SSN resolutions, but radar cannot have less than two-way 


delay or large variations regardless of processing errors. 

The excess delay equals light time for the instantaneous 
range, and the residuals match the radial distance that the 
spacecraft would travel in that time. Corresponding shifts 
in the telemetry spectra are implied by the consistency of 
the demodulated range codes with the delay in the carrier 
affecting the DSN Doppler. Both effects are traced to large 
radial Doppler rates not seen with orbiting satellites; their 
general absence beyond orbit range is also explained below 
by the spectral selection in the receiving process. 

The core contribution, with the broadest significance, is 
the explanation of the delays and the shifts themselves as 
properties of travelling wave chirp spectra, since they are 
impossible from traditional sinusoidal spectra. The chirps 
relate to sinusoidal wave spectra as rotations over the local 
frequency-time planes at the source and the receiver, the 
rotated frequency axes signifying phase accelerations, and 
equivalently clock accelerations. The shifts result due to 
causality and the finite speed of light, whose manifestation 
in the rotated view resembles expansion of space. 

The result finally reveals, and closes, a fine gap between 
d’Alembert’s general solutions and Bernoulli’s solution to 
the vibrating string problem as a series in sinusoidal waves 
[241125], that has been thought complete because of Fourier 
theory, but makes sinusoidal transport look fundamental 
and special. The constancy of frequencies is often assumed 
as sinusoidal wave solutions (cf. [26] §1.3], (27] §10-8]), or 
obtained as eigenfunctions of time invariant Hamiltonians 
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(cf. [27] §10-3,4], [251 §28,29]). However, the stationarity 
of source dynamics or constancy of carrier frequencies has 
no bearing on decomposition at a receiver, which is strictly 
computational and dictates the spectral components seen, 
and thus also lags in time varying component properties, 
including frequency and wavelength in chirps, which must 
be travel invariant to satisfy d’Alembert’s equation. 

The shifts then arise as chirp lags, but empirical proofs 
were needed for both the computational choice and reality 
of the lags, since the distance information is impossible per 
current theory. The computational aspect and availability 
of distance information in waves are specifically proved by 
the absence of the anomaly in the ESA Doppler analysis, 
which uses a Fourier transform [29] . in the Rosetta 2009 
flyby [32] 5 while its presence in range data from the same 
signal, demodulated using a carrier reconstructed with the 
Doppler rate, required a false ephemeris correction ED- 

The SSN range residuals in the 1998 NEAR flyby and 
the ephemeris discrepancy in the Rosetta 2009 flyby thus 
bear a fundamental significance complementing relativity, 
of distinguishing a spectral reference frame from physical 
space-time. The distinction decouples the wavelengths of 
reception or observation from the source spectrum, since 
the received spectrum can be arbitrarily shifted by suitable 
choice of chirp frequency rates for any source distance, so 
that tera-hertz or X-ray images can now be obtained under 
visible illumination, for example. In communication, the 
capacity of a channel is similarly defined by the sinusoidal 
assumption |32] , but signals of arbitrary wavelengths could 
be received simultaneously as chirp modes, by using shifts 
to place them in the transmission band of the same optical 
fibre, whose capacity would be then unlimited [33]f34]. 

The SSN residuals and their implication of excess delay 
in DSN and ESA data, are explained in the next section. 
The theory of chirp travelling wave spectra is given next, 
followed by quantitative analyses of the SSN residuals and 
the flyby anomaly, substantiating the above. 

Indication of the excess delay in DSN data. To 

an observer using an accelerating clock, a sinusoidal wave 
should appear as a chirp having the reverse rate of change 
of frequency, and chirps with the same frequency rate, as 
sinusoids. Chirps waves necessarily exhibit frequency lags 
that yield range in continuous wave frequency modulated 
(CW-FM) radars. In the accelerated clock view, the chirp 
lags would appear as frequency shifts which are impossible 
in sinusoidal waves, and the shifts would be proportional 
to travel, inconsistent with wave propagation as currently 
known. Fig. [T]is the graph of the SSN residuals reproduced 
from [14] , with scales of distance and one-way travel times 
for light overlaid to expose their distance proportionality. 

The 900 m residual at the start of tracking by the SSN is 
exactly the range error that would occur in about 131 ms, 
representing an optical path length of 33, 000 km, at the 
radial speed of 6.870 km s -1 , and it far exceeds the known 
resolutions of 15-25 m at Altair and 5 m at Millstone SSN 
stations. The negative sign is from the original graph and 


can only mean that the NEAR spacecraft was that much 
closer, according to SSN radars, than estimated by DSN. 

The delays are also too large to blame radar processing. 
Coherent radars perform phase correlated integration only 
to extract weak echoes over noise. The radar use of echoes 
for round trip timing eliminates ambiguities of modulated 
range codes, which get repeated and are periodic, but are 
the source of DSN and ESA range data. The SSN datasets 
thus denote true round trip times, and large errors solely 
during flybys would be in any case unlikely. Occam’s razor 
dictates, given the negative sign, that the DSN signal had 
an excess delay impossible by current ideas, but consistent 
with chirping due to acceleration, as follows. 



Fig. 1: SSN residuals from [14] with lag, range annotations 

Denoting the instantaneous range errors as Ar, and the 
radial speed as v 0 , the lag times in the figure are given by 
At = A r/v a , and the one-way ranges, by r = c At — r e « 
c Ar/v 0 — r ei where r e = 6.371 km, the earth’s radius. The 
slope of the residuals thus signifies proportionality of the 
range error to travel time as Ar = vAt. The consistency of 
the DSN Doppler and differenced range data EE33 implies 
the same error affected the DSN Doppler. 

The nonrelativistic two-way Doppler is given by Av = 
2 vvIc at frequency v for a velocity v = dr/dt , so the DSN 
phase counters yielded smaller shifts Av' = 2 vv'/c < Av. 
The observed travel time proportionality more specifically 
implies velocity error d(Ar)/dt = Av = d{vAt)/dt = a At, 
where a is the approach acceleration. A Doppler lag can be 
only significant during a Doppler rate d(Av)/dt = 2 va/c, 
whose lag d{Av)/dt x At would be therefore of frequency. 

The uplink frequency was ramped to keep the downlink 
steady during the flyby [3lj, so the delay and lags occurred 
in the uplink, and were carried into the downlink by the 
phase-synchronous transponders onboard (cf. [I] §III-A]). 

The DSN carrier loop is designed to track the downlink 
carrier frequency continuously even when its Doppler shift 
is changing (cf. [33] [lj §111]), hence the DSN phase counts 
are of cycles of changing periods, whereas Doppler theory 
was formulated for change in sinusoidal wave periods [36]. 
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The ESA extracts the Doppler using a Fourier transform 
|29j , and thereby conforms to the sinusoidal definition even 
during accelerations, since each output “bin” of a Fourier 
transform is a count of cycles around a single frequency. 
The bound of 4 pm s _1 ±44 pm s _ 1 (ler) stated against the 
anomaly in Rosetta’s 2009 flyby [3HI are just the resolution 
and phase noise in the ESA’s Fourier transform. 

The reconstructed carrier used for demodulation had to 
have been again a chirp, however, given the Doppler rate. 
As Rosetta approached the earth along its orbital motion 
from behind for gravitational boost (see [37] for all three 
flyby trajectory diagrams), the earth would have receded 
over the excess delay in the range data. The 13.34 km s -1 
perigee velocity and 2483 km altitude suggest 8 ms excess 
delay, and 110.5 m range error as the magnitude of ESA’s 
erroneous ephemeris correction. 

In CW-FM radar, the frequency lags yielding the range 
comprise cumulative change of transmitter frequency over 
the radar pulse round trips. Although the Doppler change 
was similarly continuous in both pre- and post-encounter 
tracking segments, and the modulated range codes yielded 
similar lags, the reception process represents a maximum 
integration time T shorter than a single bit in a modulated 
range code, so the implied lags and frequency rate of the 
modulation side-band spectrum, cannot have depended on 
integration through the round trip. That is, lags in a chirp 
spectrum depend only on the instantaneous rate, and not a 
cumulative change of frequency, unlike cosmological shifts. 

The residuals are thus evidence for chirp spectra bearing 
lags exceeding the total carrier variation over the receiver 
integration times T, and for realizability of fractional lags 
2 = Iw/v ss j3r/c /3T, variation of the receiver local 
oscillator (LO), which followed the Doppler rate. 

Chirp travelling wave spectra. — The general form 
of d’Alembertian solutions /(r±ct) requires / invariant of 
the retarded time (t — r/c ). Invariance in t or r separately, 
generally assumed for separating space and time parts of 
dynamical equations, would be redundant for waves as the 
d’Alembertian solutions are already most general. Rather, 
as characteristic solutions defined by and for the constraint 
of constant frequencies, sinusoidal waves were never most 
general. The assumption of constancy avoided a problem, 
however, that any variation of frequencies with distance r 
or time t would make the received waves differ from those 
observable at the source, i.e., at t = r = 0. 

Yet, any travel-invariant property £ of a travelling wave, 
hence other than amplitude or phase, should be allowed to 
vary over time locally at points on the wave path, and must 
then exhibit the lags A£ = £(f) — £(t — r/c) = £(t) — £r/c+ 
£(r/c) 2 /2! — ... = t;(t)[l—/3r/c+p( 1 \r/c) 2 /2\ — ...], where 
P = r'dS/dt, /3 (1) = t~ l d 2 £,/dt 2 , pW = £- x d 3 £/dt 3 
are fractional derivatives of £ by the receiver’s clock. This 
is unlike the Hubble shifts, which are characterized using 
proper time along the path in current theory. 

Fig.[2]shows that such lags must occur in the wavelength 
of a chirp wave because its local value around each crest 


and trough moves with the wave. The lag AA = (A 4 — Ai) 
at time £2 at receiver R must occur, in a locally measurable 
sense explained ahead, as the waveform stays unchanged 
by travel. The fractional shifts z ~ 1 — Pr/c additionally 
imply time dilations, via the Fourier inverse 

+ (i) 

the amplitude factor denoting stretching of the energy over 
a dilated interval. Equation m governs all uniform shifts, 
including both Hubble shifts and Doppler, as highlighted 
recently by the Cassini-Huygens link failure as the signal 
dilation was overlooked [35] . Dilations were not considered 
in Dirichlet’s conditions, which assured the completeness 
of Fourier theory [24][25]. As a receiver’s local oscillators 
can be independently varied at arbitrary fractional rates /?, 
and would yield the corresponding chirp spectra as proved 
ahead, the reconstructed waveforms would differ from the 
arriving waves by arbitrary time dilations, which further 
depend on the distances of the individual sources! 



Fig. 2: Lags in chirp waves 


As a prediction in a differential form from radar imaging 
mm, this had made no sense and seemed causally flawed 
[39j . It is finally explained by the computational character 
of a chirp spectrum in Fig. [3] as a rotation of the receiver’s 
local frequency ( RRq ) and time (— RRt ) axes, denoting 
the local evolution of the spectral components in time by 
the receiver’s clock. The constant frequency of a sinusoidal 
component would be represented by vertical lines like BC. 
The inclined lines GC , HF denote chirp components with 
frequencies increasing over time. The spectrum at present 
time t 2 is represented by the same coefficient values on the 
frequency axis RRq regardless of the inclination. 

With the inclination, however, excess one-way delays are 
incurred, just as in the DSN Doppler, that result in shifts 
exactly equal to cumulative change from an earlier state at 
the source, so the distance information bears the penalty 
of excess delay. Each chirp line, projected indefinitely, not 
only attains every possible frequency at some instant, but 
is identical to every other chirp of the same inclination by a 
simple displacement in time. This equivalence leads to the 
excess delay, as the travel delay acts against the frequency 
change. Conversely, were the angle of inclination LCBF = 
tan _ 1 (|CF|/|HC'|) = tan - 1 (/3At/Af) = tan -1 /? made 0, 
the chirps would become degenerate vertical lines through 
C and F that overlap no longer if displaced in time, so the 
delay and the distance information both vanish. 

These details, and relations to causality and the speed 
of light, are revealed by incorporating the source frequency 
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(SSq) and time (— SSt ) axes, with corresponding source 
chirp lines JD and AE parallel to GC and HF. Sinusoidal 
transport would be represented by parallel lines like DC 
and EF connecting equal values on the source and receiver 
frequency axes. Hubble’s law would require inclined lines 
like EC to produce shifts Aw = W4 —Wi = \CF\ at distance 
r and Awi = \LM\ ss \CF\ri/r at distance rq = \EM\. 



Chirp spectra would be thus time invariant like Fourier 
spectra, but exhibit distance proportional shift factors and 
dilations with the receiver’s choice of j3 and its derivatives, 
because the chirp spectra start fully shifted and dilated at 
the source! The total energy is also clearly unchanged. 

The inclined axis RpRp,n, denotes the chirped spectral 
view, given by the DSN and ESA range data during flybys, 
in which local chirp histories GC and HF seem normal to 
the frequency axis, but travel lines AC, KL unaccountably 
seem inclined. The inclination of axis is equivalent to the 
receiver’s clock acceleration inclining the components; the 
segment | PF denotes the relative phase accelerations \PF\ 
that are not apparent in the rotated “reference frame”, in 
which the chirps appear as a Fourier spectrum with shifts 
|CP| ss Aw, due to skewing of all travel lines AC, KL to 
longer wavelengths, as if space itself were expanding. 

Reception and orthogonality. — In any frequency 
modulation scheme, including phase shift keying (PSK) in 
deep space telemetry m, can be described by a random 
variable Ll m denoting the instantaneous modulation. Both 
at the DSN receiver and the spacecraft transponder, the 
carrier loop phase locks imply, upon allowing for frequency 
variations, the first order product integral condition 


The r—ct invariance required of d’Alembertian solutions 
more particularly requires lines like AC and KL inclined 
at LDAC = tan _ 1 (|DC|/|AD|) = tan _ 1 (r/Af) = tan -1 c 
with respect to the time axes, but parallel to the distance 
vector SR, so as to connect equal component frequencies of 
the source current and receiver voltage spectra, regardless 
of whether the connected frequencies belong to chirps, as 
denoted by lines IE and HF, or to sinusoids, represented 
by lines AD and BC, respectively. The inclination denotes 
wave speeds c < 00 , and is along of increasing time from 
source ( A ) to receiver (C), conforming to causality. 

More importantly, a component with angular frequency 
wi at C on chirp line GC should correspond to the same 
angular frequency wi in source history (A), but belong on 
chirp line IE that changes to W 4 by time £2 (E). However, 
an atom emitting at angular frequency W 4 at £2 ( E ) would 
have been observed locally at W 4 also at t 0 ( N ), and the 
same should hold for a steady carrier transmission. It thus 
appears that the d’Alembertian travel lines like AC either 
require amplitudes to shift with travel, from N to C, which 
would conflict with the d’Alembertian invariance; or chirp 
spectral decompositions, which can only produce inclined 
histories like HF and AE, must be impossible, so the lags 
Aw would require nonlocal simultaneous measurements at 
source and receiver at £ 2 - The second case is cannot hold 
since the inclinations /3 could be infinitesimally small, and 
the decomposition is in any case purely computational. 

The answer is that the construction already implies that 
at nonzero j3, wi is seen only at distances r = c Aw//3. The 
amplitude at C comes from A, which is precursor to E at to 
and to N at to - At. The chirp spectrum thus reconstructs 
distributions at the past times to — At = t 2 — 2At, where 
the factor 2 relates to the excess delay. 



^io c + n m )e 0 '^- r ^ 


exp 
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dt 


— 27 tS (w c + (D m ) - ui 0 )6(P' - /3), 
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where w c is the carrier; w Q is the loop voltage-controlled 
oscillator (VCO) frequency; p and f3' are fractional rates of 
the VCO and a received spectral component, respectively; 
and T is the loop filter time constant. T is set below 1 Hz 
in DSN carrier loops in order to suppress both phase noise 
and modulation [35]. The /3 _1 factor is from integrating 
the exponential chirp w(t) = woe^* for get the phase, and 
vanishes in the phase derivative via L’Hopital’s rule. 

Equation © constitutes the orthogonality condition for 
exponential chirp waves without modulation (fl m = 0 ), 
and including the case of T —> 00 , since a travelling wave 
of the same instantaneous frequency and rate of change as 
the receiver’s LO (~ w D ) contributes in every cycle to the 
integration performed by subsequent filters, but any other 
component contributes over at most a cycle. In a Fourier 
transform, nonmatching components contribute at every 
few cycles indefinitely, so Fourier convergence depends on 
Cesaro means, and is weaker in this sense. 

The orthogonality looks weak for distinguishing between 
say, a chirp of fractional rate 1 s _1 at 5 THz from a 5 THz 
sinusoid as their phases would differ by only 10 - ' rad over 
10 5 cycles, but in a spectral selection or decomposition, all 
families of curves over local frequency-time planes (Fig. [3]) 
must be assumed available. The contributions from /3F5/3 
pairs then cancel out for 5/3 7 ^ 0 , just as in the interference 
of alternative paths in Fermat’s principle. For decoding or 
demodulation, equation © relates the modulated carrier 
and LO statistically over shorter integration times T for 
the modulation bandwidth, assuming (S 2 m ) = 0 , since the 
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d.c. (direct current) is suppressed in deep space telemetry. 
The consistency of the DSN range data with its Doppler 
implies that the modulated chirp spectrum had the correct 
phase offsets fl m relative to the chirp carrier. 

Larger lags, of 9 MHz at 1 AU for the same acceleration 
0.5 m s -2 as at loss of signal (LOS), would shift the chirps 
out of the filter pass-bands, so the signal presumably gets 
demodulated from the Fourier spectrum without lag^j. 

Explanation of the SSN residuals. — The net gain 
in speed was only (6.87 — 6.83)/6.87 « 0.6% 0 Fig. 3a], 
with most of the acceleration close to earth after the SSN 
tracking in increasingly tangential motion. The uniformity 
of the 10 min ticks in the equatorial view [9] Fig. 1] and of 
similar ticks in the north polar view 1141 Fig. 9], which are 
expanded due to projection, suggest that the mean speed 
v a = 6.85 km s -1 would be adequate for present purposes. 

The 219 min gap in the DSN tracking then represents 
6.851 km s -1 x 219 min = 90, 000 km of trajectory. LOS 
occurred lh 8 min before periapsis and acquisition of signal 
(AOS), at Canberra, at 2h 31min after periapsis, so the 
range was 90, 000 km x 68/219 ~ 27,950 km at LOS and 
62, 070 km at AOS. Tracking at Altair ended 36 min past 
LOS at 06:51:08 and had started at 06:14:28, for a total of 
2200 s, so the tracking started 4120 s before periapsis, at 
r « 4120 x 6.851 km s -1 = 28, 226 km. The one-way delay 
was therefore At = —r/c « —94 ms, implying range error 
e r = v At = —94 ms x 6.851 km s -1 « —645 m, about 25% 
smaller than in Fig. [T] The error decreased with the range 
rate at de r /dt ps vd(At)/dt = v 2 /c = (6.851 km s _1 ) 2 /c ~ 
0.313 m s -1 , over 1187 s from 06:25:25 to 06:45:12, hence 
by 0.313 x 1187 « 186 m, which is within 10% of the 200 m 
Millstone decrease. Fig.[l]shows two sets of residuals, since 
they are projections of the same lag of the (pre-LOS based) 
DSN estimate behind the true trajectory in the direction 
of each of the two SSN stations. The ground track diagram 
p~4l Fig. 7] shows the trajectory pointed towards Millstone 
initially, implying a faster initial decrease of range, hence 
greater initial values for Millstone, as seen in Fig. [T] 

Explanation of the flyby anomaly. The delay 
means that DSN underestimates pre-encounter approach 
speed and overestimates post-encounter recession and thus 
infers an anomalous velocity gain in earth flybys whenever 
the tracking is discontinuous across periapsis. If tracked 
continuously, however, the delay in the Doppler’s change 
of sign at periapsis (Fig. U) should cause a negative Av. 

The negative Av in Galileo’s second flyby was concluded 
from around periapsis, since it was at first thought masked 
by atmospheric drag HUH- The tracking was unbroken 
in Cassini’s flyby that also showed negative Av m- 

As the excess delay varies with range, the true velocity 

1 The lags should be about 18 Hz at lunar range, but inband 
chirps would then face echo suppression due to their 1.2 s excess 
delay. A carrier loop lock to the Fourier spectrum would produce a 
piecewise frequency approximation of the Doppler rate as the VCO 
carrier, so each range code bit then gets retrieved from the Fourier 
spectrum. 


profile, given by the differenced SSN range, and the DSN 
Doppler would be closest at periapsis, and cannot really be 
parallel. The slopes of the residuals were thus “irreducible 
through velocity estimation” m, though both curves were 
monotonic over the SSN tracking period, as shown. 


6.87 km/s 


- - DSN Doppler 
true velocity 



^ Av < 0 6 ' 83 km/s 

Fig. 4: DSN Doppler and its lags during flybys 


The velocity error at AOS should cause post-encounter 
data to be inconsistent with the pre-encounter trajectory, 
and vice versa. Acceleration due to earth’s gravity at AOS 
range would be a = 0.103 m s~ 2 , implying a velocity error 
Av = —aAt ps 21.4 mm s _1 , or a 603 mHz ~ 10.7 mm s -1 
Doppler amplitude at the downlink frequency. These are 
about 20% of the reported 760 mHz = 13.5 mm s -1 [9]. 

Canberra’s latitude of 35.2828° means it is 6371 km x 
cos(35.2828°) = 5201 km off the earth’s axis. The —71.96° 
declination of the post-encounter velocity asymptote then 
implies 5201 km x cos(71.96°) « 1611 km of diurnal range 
and (1611/62070) x 603 mHz « 15.6 mHz diurnal Doppler 
oscillations. The larger actual 50 mHz amplitude is due 
to a smaller declination at AOS, and to a misprediction of 
direction [9], possibly worsened by the error at LO^|. 

The velocity error Av = —aAt = —ar/c also explains 
the r _1 decay in the post-encounter oscillation graphs [9] 
H, since oocr 2 due to earth’s gravity. 

Conclusion. — All features of the flyby anomaly are 
thus explained by a delay proportional to range in Doppler 
and range data derived from the telemetry signal that was 
chirped due to the Doppler rate, and should be impossible 
by current notions of wave propagation. 

More particularly, the distance proportionality of the 
delay and of the equivalent frequency lags in the telemetry 
spectrum are given by two independent radar tracks which 
had been overlooked for over a decade, oddly, in the very 
quest for deviations from the known laws. 

The chirping and lags should be also easy to verify over 
ground distances at radio frequencies, with no motion or 
the difficulties of optical implementation (cf. j34l). 
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